o 

(N 



> 

o 



o 



X 



Nearly Doubling the Throughput of Multiuser 

MIMO Systems Using Codebook Tailored 

Limited Feedback Protocol 

G. Wunder, J. Schreck, P. Jung 
Fraunhofer Heinrich Hertz Institute, Berlin 



Abstract 

We present a new robust feedback and transmit strategy for multiuser MIMO downlink commu- 
tyj ' nication systems, termed Rate Approximation (RA), and analyze its performance. The new scheme 

combines flexibility and robustness needed for reliable communications with the user terminal under a 



limited feedback constraint. It responds to two important observations: One is that it is not so significant 
to approximate the channel but rather the (potential) rate itself so as to mimic the optimal scheduling 
decision at the base station. The second observation is that a fixed transmit codebook at the transmitter 
^N I is often better when simultaneously the channel state information is more accurate. Both observations 



are incorporated in the new scheme where the transmit and feedback codebook are strictly separated 
and user rates are delivered to the base station subject to a controlled uniform error regardless what 
the scheduling decision is. The scheme is analyzed and proved to have better performance below a 
certain interference plus noise margin and better scaling properties than the classical Jindal formula 
ii , when considered in the very same setting. LTE system simulations sustain the analytic results showing 

performance gains of up to 50% compared to zeroforcing when using multiple antennas at both the 
base station and the terminal, and up to 70% when using single antennas at the terminal. Finally, a new 
feedback protocol is developed which inherently considers the transmit codebook and which is able to 
deal with the complexity issue at the terminal. 

I. Introduction 

Multiuser multiple input multiple output (MU-MIMO) communication systems have been 
in the focus of intensive research over many years. The optimal transmission technique for 
these systems is dirty paper coding (DPC) |[I1 which, under perfect channel state information at 
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the transmitter (CSIT), achieves superior performance gains over linear schemes. However, in 
practical systems CSIT is obtained via a rate-constrained feedback channel, which is known to 
be a sensitive part of the overall system and must be carefully designed. An extensive overview 
of so-called limited feedback MU-MIMO systems can be found in [2J. References ^, JH set 
the standard for performance evaluation in MU-MIMO systems. 

In this paper we revisit the limited feedback problem in MU-MIMO systems. We consider 
linear beamforming techniques and assume that the transmit beamforming vectors are defined 
by a fixed transmit codebook. In contrast to previous work we use a different codebook for 
the feedback and apply a new feedback strategy which we call Rate Approximation (RA). The 
main idea is that the terminal selects a channel quantization vector from a feedback codebook 
considering any possible combination of beamforming vectors from the transmit codebook. As 
we will show, this enables the base station to approximate the user rates subject to a small 
uniform a priori error. Given the feedback message, the base station is then permitted to assert 
any beamforming vector from the transmit codebook for some optimization purpose (not just 
the beamforming vector dictated by the user if scheduled). 

Let us provide a striking example: Consider a setup where the users are served on unitary 
beamforming vectors. It is well known that such a scheme achieves the optimal throughput for a 
large number of users ^. However, it has been mostly overlooked yet that the scheme performs 
excellent for a finite number of users as well, in particular with limited feedback under the new 
RA strategy. This is illustrated in Figure [IJ In the analytical section (see Corollary [I]) we make 
exact this observations precise in this paper. 

Organization and Main Results: In Section |II] we introduce the system model and in Section 
Uni the RA scheme is introduced. In Section |IV] we analyze the performance of the RA scheme. 
Our baseline is ZF and the operating point is such that the base station always serves as many 
(independent) users as transmit antennas. This model is essentially in accordance with f3] and 
enables us to make stringent comparisons. The analysis is performed in three steps: 

1) We consider perfect CSIT and show that unitary beamforming (UB) can achieve a signifi- 
cant performance gain over ZF for a large fraction of a practically relevant signal-to-noise 
ratio (SNR) range. 

2) We analyze the a priori rate error at the base station (i.e. before any scheduling decision) 
for each individual terminal evoked by our RA feedback strategy and prove that it is not 
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Fig. 1. Spectral efficency vs. average SNR: comparison of ZF with perfect CSIT and beamforming under a fixed codebook 
with partial CSIT due to Rate Approximation. The bounds for ZF are given by (3j. It is observed that RA not only performs 
better for realistic SNR but also is immensely robust to channel uncertainty. 



only typically much smaller than that of ZF but has also better scaling properties compared 

to the classical result by [I3]| which is even improving in the number of transmit antennas. 

Consequently, together with 1) we characterize the single cell SNR operating points where 

a fixed codebook has better performance than ZF, clearly setting a new standard. 

3) We outline an advanced vector quantization problem for the RA scheme replacing the 

common chordal distance with a new distance function which inherently uses the structure 

of the transmit codebook. 

In Section |V] we underline our results with LTE system simulations showing the benefit 

obtained by the proposed RA scheme and develop a suboptimal feedback protocol dealing with 

the complexity issue. This feedback protocol is proposed to replace the common approach for 
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LTE. Finally, in Sec. |VI]the conclusion is drawn with emphasis on the impact on future standards. 
Notation: Bold letters denote column vectors and bold capital letters matrices. The inner 
product between vectors a and b is defined as (a, b) = a^b, where a^ is the conjugate transpose 
of the vector a. §""^ is the unit sphere in C" for < n G N; a^ ^ denotes the j-th component 
of the vector a^ and [a.i]j^i „ = [o-i, 0,2, . . . , a„] is the matrix with column vectors a^. The 
non-negative integers are denoted as N+. 

IL System Setup 

A. System Model 

We consider the MU-MIMO downlink channel of a cellular system where a base station, 
equipped with Ut transmit antennas, serves multiple users, equipped with Ur receive antennas, 
on the same time and frequency resource with a single data stream. The users are collected in the 
set U. Let x G C"'^^ be the transmitted signal. User m receives the transmitted signal through 
the channel H^ G C""^^"* and applies a receive filter u^ G C"''^^ to recover its intended signal, 

Vm \""m) -l^m'^/ ~r ^m \''"mi '^/ ~r ^m 

where n^ = {um, rim) ~ CJ\f{0, a^) is additive white Gaussian noise. The effective channel 
from the base station to user m is given by h„^ = {um, Hm)\ the normalized effective channel 
will be denoted by hm = /im/ll^mlh and the effective channel norm is /i^ = ll^mlh- Further, 
we assume each user m has perfect knowledge of its own channel Hm and that the channels 
are constant over one transmission interval; no fading model is imposed. 

In each transmission interval the base station selects a subset S CU of users for transmission 
and assigns each user m E S a beamforming vector out of a transmit codebook C = {wi, W2, ...}, 
known to the base station and all users a priori. Sometimes, we will also assume that the codebook 
elements constitute an orthonormal basis (ONB) of C"', such that for every / G C"* 

rtt 

In that case we speak of UB. The assignment of users to beamforming vectors is defined by a 
mapping 

vr : 5 ^ C, 
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that maps each element in 5 to a codebook element from C. We assume that |iS| < Ug < rit, 
where Ug is the maximum number of scheduled users per resource. In the sequel we do not state 
the domain of vr explicitly, if it is clear from the context. Define the complex information symbols 
intended for user m as dm. € C, then the transmitted signal can be given by the superposition 




In the following the power allocation is uniform, that is, the base station distributes its available 
power equally among all users. The transmitted signal must satisfy an average sum power 
constraint, i.e. E [||a;||2] < P. Define 



,2 _ -* A'-rra 



as the receive SNR of user m eU, which is a measurable quantity at the receiver. The achieved 
sum rate under any user selection S and any mapping tt : 5 — )■ C is 

R (tt, H) = ^ rmiir, Xmhm) (2) 

where the effective channels are defined by the composite matrix H = [\mhm]mes for some 
set of users S CU. The per user contributions to the sum rate are given by the Shannon rates 

r„(7r, Xmhm) = log M + -^^ rnm,^{m.) _ \ ^^^ 

where the effective channel gains are defined as 

Throughout the paper we assume maximum sum rate scheduling, for instance, with perfect CSIT 
the optimal scheduling decision tt// : 5 — )■ C is 

[S, TXu) = arg niax R (vr, H) . (4) 

scu 

However, due to the rate-constrained feedback channel, the base station takes its decisions based 
solely on partial CSIT from each user. Clearly, these decisions should match the optimal decision 
as good as possible. This challenge can be most efficiently approached by the following RA 
scheme. 
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III. Rate Approximation 

Before introducing and analyzing the RA scheme we review some basic results for ZF, which 
is our baseline performance. 

A. Baseline Performance: ZF 

For the ZF scheme with perfect CSIT, each beamforming vector WzF,m £ C"' is chosen to 
be in the null space of the matrix [hi]i^gxr^x (see e.g. Q). The sum rate is then given by 

RzF (H) = Y, log (1 + \l<Pt,ZF) , (5) 

r?ig5 

which is almost surely non-zero for any receive SNR Am > 0. Here, we have defined 0^ ^f = 
\{hm, 'WzF,m)\'^- Furthermore, define R{H, Vrvq) as the sum rate of ZF under limited feedback 
based on random vector quantization (RVQ) [16], where Vrvq = {i^^^K . . . , i/'^ '} is a codebook 
of random vectors isotropically distributed on the unit sphere. In contrast to the RA scheme 
(introduced later on) user m chooses his feedback message to minimize the chordal distance 



dc{h^,u) = ^/l-\{i^,h,,)\^. (6) 

From [|3l Theorem 2] we have the following performance bound for ZF with RVQ. Suppose 
that hrn,i ~ CJ\f{0, 1) independent across users and antennas. Then, limited feedback with B 
feedback bits per user incurs a throughput loss relative to ZF with perfect CSIT according to 

Ai?RVQ = ^H [RzF (H)] - Eny [R (H, Vrvq)] 

P B_ 

<nilog(l + ^2 -t-i). (7) 

B. Rate Approximation Scheme 

Assume the base station has partial CSIT (^m, Vm) for each user m eU. The channel direction 
information (CDI) !/,„ G V is given by an element of a feedback codebook V = {v>^^\ . . . , v^'^ ^} 
that consists of a collection of normalized vectors i/W G S"*^^ and is a priori known to all users 
and the base station. The scalar ^^ G M can be interpreted as the channel quality information 
(CQI). In the sequel we assume that the CQI is perfectly transferred to the base station, which is 
a typical assumption. A natural choice for the CQI is 'dm = Am if CDI can be perfectly tracked, 
but under partial CSIT the choice is in general crucial. 
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If the beamforming vectors are restricted to a fixed codebook C tlie sclieduling decision based 
on partial CSIT V = [Omi^mlmes '^^^ ^^ found by solving 

(S, TTy) = arg max R (vr, V) , (8) 



scu 



which is a combinatorial problem that can be solved either by a brute force search over the user 
sets SCU, with |iS| < n^ and the mappings vr : 5 — !■ C or, alternatively, more efficiently in a 
greedy fashion DTI, flSll. 

Let us now derive the RA scheme. The aim of the RA scheme is that the scheduling decision 
under partial CSIT matches with the optimal scheduling decision (Hj) as good as possible. 
In the following we demonstrate how the error can be calculated thereby circumventing the 
optimal combinatorial solution. Define the average rate gap between ZF with perfect CSIT and 
beamforming based on a fixed codebook with perfect CSIT as 

Ai?csiT := Eh [Rzf (H)] - Eh [R {tth, H)] 

and the average rate gap between the real sum rates R (vr, H) and the approximated sum rates 

R (vr, V) as 

AR {it) := Eh [R (tt, if)] - Eh [R (tt, V)] . (9) 

Now, the rate gap between ZF with perfect CSIT and beamforming based on a fixed codebook 
with partial CSIT can be bounded from above by 

ARzF = Eh [Rzf (H)] - Eh [R (ttv, H)] 

= A/^csiT + Eh [R (vth, H)] - Eh [R (vry, H)] 

= Ai?csiT + AR (tth) + Eh [R (vr^, V)] - Eh [R (vry, H)] 

< Ai?csiT + AR (tt^) - ARinv) (10) 



< AiicsiT + 2E 



H 



Emax |r„(7r, XraKa) - r^(vr, ?9„z/, 
Om, Co 



*5j7],CtS 



(11) 



where (iTOl) must hold since vry is the optimal mapping of users to beamforming vectors under 
channels V. In (fTTI) we defined the set 

Sm ■■= {{m}, {m, m + 1}, . . . , {m, m + 1, . . . , m + n^ - 1}} , 
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that is, the set of possible user selections which include user m. Moreover, we exploited that 
the rate gap AR (tih) — AR{ttv) is bounded from above by the worst case rate gap 



ARra := 2 ■ E 



RA •= -i ■ li^H 



Om Co 






(12) 



From (fTTI) we observe that to bound ARzp it is sufficient if each user individually minimizes 
its individual rate gaps |rm(vr, A^^m) — ^ml^i", ''^m^m)| for any Sm- This observation is the 
fundament for the following RA scheme. 

To determine its feedback message each user m E U must find a tuple {'dm, t^m) G (^, V) 
that minimizes the RA distance ' 



d{Xmhm,^mi^m) = Hiax Ir^ (tt, A^/i^) - r„ (vT, ^^I/m) | . (13) 

That is, each user m eU finds its feedback message by solving 

{^m,l^m) = SiTgmmd{Xmhm,^J^) ■ (14) 






The RA scheme can be easily extended to users with multiple receive antennas n^ > 1. In this 
case for each scheduling decision n : Sm -> C the optimal receive filter can be considered in 
the RA distance, i.e. the real rates are given by 

Tmi^r, Kn) = max Tm U, {u, H m)") ■ 

Although not apparent at this point let us formulate some decisive advantages of the RA scheme: 
first, in the RA distance c?(-, ■) the transmit codebook matters which seems good engineering 
practice as we use all the available information. Second, the terminals provide an uniform error 
which indicates how well the rates are approximated and leads to inherent robustness. This 
becomes particularly beneficial in the LTE multi antenna case where CSI information is averaged 
over the subcarriers (see Simulations in Section JV]). Third, the RA scheme is amendable to 
codebook optimization [|9l due to our new distance function. Finally, let us establish that the RA 
distance is indeed different compared to the chordal distance by the following example. 

Toy Example. To illustrate how the RA scheme operates we consider a toy example in M^. 
We compare the feedback decisions taken under the RA distance (fT3l) and minimum chordal 



A closer look reveals that it is neither in all cases a distance on C"' nor on the Grassmann manifold. 
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b) 



Fig. 2. Sets TZi and Mi in the toy example, a) RA distance: If a user experiences a channel h G TZi the quantized channel 
Vi is fed back, b) Minimum chordal distance: A user experiencing a chaimel h £ Mi feeds back channel quantization Vi. 



distance ^ by user m E lA. Assume ng = 2, nt = 2 and rir = 1 such that the channel 
vector is given by hm G IR^- The transmit codebook is given by the columns of the identity 
matrix C = {(1 0)^, (0 1)"^} and the feedback codebook is given by a rotated version of the 
transmit codebook. The CQI is equal to the receive SNR ^^ = A,^. For the RA distance we 
define the sets 

TZi= IhmeR^ : i/W = argminrf(A„^„, A^i/) 

That is, all channels that result in the codebook element Vi E V. In a similar manner for the 
minimum chordal distance we define the sets 



Mi = <hjn ER : i^-'*' = a.Tgmm dc{hm, i^] 

Figure |2] a) shows the two sets IZi and 7^2 and Figure |2] b) shows the sets M.i and M.2- 
We observe that the RA distance results in a different CDI feedback decision compared to the 
minimum chordal distance. The RA feedback decision is obvious more oriented on the transmit 
codebook. 

In the next section we show that Ai?RA can be made small and that there is an SNR range 
for which Ai?csiT is negative. 
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IV. Performance Analysis 

A. Problem Formulation 

Many papers prove that a particular transmission scheme achieves the optimal multiuser 
multiplexing gain. That is, for sufficiently large \U\ the sum rate scales like nt log log |W|. For 
instance this was shown for random beamforming (Si, ZF [[TOll and UB IfTTI . However, since 
rates and the number of users are finite in a practical system, the significance of these asymptotic 
results can at least be questioned. Putting it the other way around: two methods achieving the 
optimal gain might behave completely different in a practical system. This implies that it does 
not say much about individual rates of users. 

Our analysis is different and more inspired by the finite user results in |[3l, [[T2l|. For the sake of 
analytical tractability, we first constrain the number of users and selected users by \U\ = \S\ = rit, 
thereby circumventing the user selection process and enabling stringent comparison to Jindal's 
results in [[3|, and also resort to UB. Later on, we will abandon this assumption and incorporate 
user selection and general codebooks as well. 

Our goal hereafter is to upperbound the rate loss with respect to ZF with perfect CSI; this is 
done in two steps: 1.) first the base station is assumed to have ideal CSIT, but operates solely on 
the set of beamforming vectors of the transmit codebook (this corresponds to evaluating Ai^csix) 
and 2.) the individual rate error due to the rate-constrained feedback channel is bounded from 
above by involving the RA scheme (this corresponds to evaluating A/^ra). 

B. ZF vs. UB: Perfect CSIT 

So far we have not used any fading model. In the following theorem, the channel is modeled 
as isotropic fading [[T3l . In particular, we fix the channel gains Hm for m E U and let the only 
randomness be an independent phase ambiguity in the channel coefficients, i.e. hm = l^mhm 
where hm E §"*^^ is an isotropically distributed random complex unit vector. Note that this 
assumption has no impact on the generality of the result; in fact it only allows us to streamline 
the presentation throughout the proof. 

The following theorem is qualitatively known but we make it mathematically more precise. 

Theorem 1. If U = S = {1,2, ... , rit}, then for isotropic fading with any non-random fii > 
/i2 > ... > fint (respectively non-random Xl> \l> ... > X^J the rate gap between ZF and UB 
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with perfect CSI is bounded from above by 

n -I- X2 \ 1 f^ I A^ ^ c(e)(l-e)log(nt-m+l) 

y^ + ^m.) t:;^- [^ + 7^7^; iTw^t 



^RcsiT < ^ log 1 + min A^ 
w/iere c (e) := (1 - e-("*-"*+i)' - ^). 

The proof can be found in Appendix lAl Theorem [T] states that for some A^ < A* the rate gap 
Ai?csiT is indeed negative; this can be seen by observing that at low SNR (i.e. A^ — )■ 0) for all 
m < rit the difference is roughly 

1 log(nt — m + 1) 



Ut-l Tit 



<0 



by choosing e.g. c (e) = ia°fn) ^^'^ ^^ large enough. Therefore, for low SNR unitary 
beamforming can indeed outperform ZF. 

C. Uniform RA Error with UB 

Let us first provide a general expression for the maximum in (fT2)) which gives us a hint how 
the RA scheme operates. Note that from now on when the RA scheme operates on a unitary 
transmit codebook we will denote this scheme by RA-UB. 

Lemma 1. IfU = S = {1,2,. ..,77,4} and C C V then for some pair h^, i^m under the RA-UB 
scheme 

inaX \rm{^T,\mhm)-rm{'^,'^r,^Vm)\ (15) 

- ^^.^°^ y + 1 + A^ (1 - max{|(/i^,«;)|M(j/^, w)\^}) 

where we defined w* by \{hm,w)\'^ < \{hm,w*)\^ for all w E C. The strategy is to pick u^ 
close to hm as long as |(i/m, tf*)P > \{hrn,w*)\'^. 

The proof can be found in Appendix |Bl While the error term is still not easily accessible the 
former lemma allows us to devise the following corollaries. 

Corollary 1. If C = V, then under the RA-UB scheme 

AR^ < f; E^ [log (l + max t-^^I^^t^^^^) 
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with CDI Um = w* and CQI defined by 

-■ i+\ui-\{h^,v^)\^r ^ 

Proof: Since, C = V we have that |(i/„,if;)p = for all w ^ w* . Plugging ([T5]) in (fT2l) 
proves the claim. D 

Corollary [T] gives a universal lower bound on the performance of any RA-UB scheme. Observe 
that a CQI similar to ( fT6l ) was proposed in (Si in the context of ZF and also in Q using random 
beamforming. Equation (fT6l) is also in accordance with [[TOl . where it was shown that the CQI 
should reflect the channel magnitude and the quality of the channel quantization. For the RA 
scheme, in the simulations, we will also consider the CQI 

dl = \l\{hr„.,u^)\\ (17) 

which can be interpreted as the effective channel of user m over the quantized channel Vm- 
Equation (fTTI) captures two important aspects. On the one hand, if the CDI is equal to the 
channel direction, the user gets no penalty (i.e. | (/im, Um)\'^ = 1) on the other hand if the CDI is 
orthogonal to the channel direction, the effective channel is zero (i.e. |(^m, ^'m)P = 0). Hence, 
the CQI (flTl) reflects the receive SNR and the quantization error. 

Another consequence of Lemma [U is that the sum rate R^ttv, H) achieves the optimal 
throughout scaling for asymptotically many users. We use the (standard) Rayleigh fading channel 
model, i.e. hm,i ~ CA/'(0, 1) for i = 1, . . . , rij and m eU. 

Corollary 2. Let the unitary transmit codebook be a subset of the feedback codebook C C V 
and each channel is Rayleigh, then under the RA-UB scheme we have 

\-^ — ^ ^ 1, \U\ -^ oo. 

nt loglog\U\ 

Proof: From |l5l Lemma 1] we know that under perfect CSIT the sum rate R{'Kh^H) 

defined in (|2l) scales as niloglog|W|. From Lemma [H we have that the rate gap ARra = 

R{nH, H) — R{'n'v, H) remains bounded when \V(\ increases. This concludes the proof. D 

So far we are still not able to effectively upperbound Ai?RA which is now settled. The following 

lemma shows that Ai?RA remains bounded when the SNR increases and that the rate error 
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depends solely on the function 



D„(B) 



where we defined Ar 



E 



H 



mm max 

i>i5m>o wec 



\m\(h,w) 



1^rn\{l^,w) 



Xr 



(18) 



At, 



- i9 



i?„ 



and B is the number of feedback bits, i.e. 2^ is the 



number of elements in the feedback codebook V. 

Lemma 2. If U = S = {1,2, ... , rif} and C C V then under the RA-UB scheme, 



^Rra < 2 5^ log 1 



m=l 



mm — 

e>0 1 



-Dn,{B) 



The proof can be found in Appendix O 

The following lemma gives a fundamental bound on Dm{B). 

Lemma 2>. If \V\ = 2^ and 'dm = ^m for all m ^ U and the transmit codebook C is unitary, 
then 

Drn{B)<c{nt)¥.H[\r]'^' 



B 

tif — l 



with 



c[nt) 



e(s^*^ 



,^j2nt-2\V{l 
1 



nt 



nt-l^ 



{fit - ly.TT '2 



C»t-1) 



and B > ^"'^ ' log[(ni — l)^/r^^^T]. For rit — l small tight bounds are known for the covering 
density BiB^'"^), e.g. e{Bj) < 1.2091 (Kershner, 1939), 0(6^) < 1.4635 (Bambah, 1954), 
Q{Bl) < 1.7655 (Delone & Ryshkov, 1963). For nt - I > 3 the Rogers bound fUB Q{Bl) < 
4(ni — 1) log(n( — 1) can be used. 

The complete proof can be found in the Appendix |Dl Note that c{nt) is close to unity but 
falls below unity not before n^ > 8 as required for improved scaling compared to Jindal's result. 
This is an artefact of the proof technique as the following illustration for the case Ut = 3 shows. 

Without loss of generality, we assume the unitary transmit codebook is given by the standard 
ONB and define the real positive vectors 

0m = {4>m,7T(l)^4>m,n{2)^^m,TT{3)) ^^^ 4^m = {4>m,7T{l)y 'i^m,n{2)^ 4>m,n{3)) 5 

per definition, this vectors have unit £i-norm, i.e., ||<^^||i = ||0mlli — ^ ^^^ define points on 



the standard 2-simplex. Further, max^ 



m,n{m) Tm,n{m) 



10^ — 4>m\\oo defines a distance 
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Fig. 3. The standard 2-simplex in 3 dimensions. The projection of the quantization points Q on the coordinate axes implies 
a worst case quantization error S. 



between two points on the standard 2-simplex. Hence, for a given feedback codebook V we can 
define the Voronoi region around the point (p'^ with u eV as, 

y(0 = {XE M3+ : \\x - cj^Uo. < \\x - 0IIIOO, V^ e V, ^ ^ u}. 

If B E {1,2,4, .. .} and rit = 3, the feedback codebook can be chosen such that the Voronoi 
regions are 2-simplices with edge length 5 < \pl. Now, using the symmetry of the covering 
and projecting the quantization points back on the coordinate axes (see Figure [3]) we get 
maxa.gy(^- ) II a? - </>,''„ ||oo = V^S = b. 

Now we can compute the volumes of the 2-simplices (the standard simplex and the scaled 
simplex) and proceed as in the proof of Lemma [3] to obtain the result 



b= max ||ic-(/)^||oo 



2 "t 
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Hence, if rij = 3 and each channel is Rayleigh, using the RA-UB scheme the rate loss due to 
the rate-constrained feedback channel scales like 

lufi I 1 -|- iiiin ■ 

m=l \ 



Ai?RA < 2Eh 



Ei°g 1 



l-A™ nt-1 



2 ^^T^' 



m=l ^ -^ 



Therefore, we have an improvement of Ut — 1 bits in the exponential term compared to Jindal's 
result for ZF with feedback based on minimizing the chordal distance (see Section ITlI-AI) . under 
the very same assumptions. 

D. User Selection and General Codebooks 

In this subsection we no longer assume unitary transmit codebooks and allow user selection 
at the base station. 

Theorem 2. If the base station selects a subset S CU = {1,2, ... ,nt} of users per transmission 
and each channel is Rayleigh, then the RA scheme with B feedback bits per user incurs a 
throughput loss relative to perfect CSIT bounded from above by 



Pm 
ARra < 4n, log ( 1 + ^Eh 



a^ 



minmax I \(h, w)\' 



U/, W) 



The proof can be found in Appendix El The expected value 

Eji" lamiaaxlKh, w)\'^ — \{i/,w)y 
uev wec 

has been shown to be analytically traceable, in the previous section, for unitary transmit code- 
books. However, its examination for arbitrary transmit codebooks C goes beyond the scope of 
this paper. It is expected that the scaling advantages of the RA scheme diminish with increasing 
transmit codebook size. But, clearly for unitary codebooks the previous result still holds which 
is remarkable, as all the 2"* possible user rates are uniformly recovered at the base station with 
better scaling properties than the classical result. 

The RA scheme is now applied in a practical scenario. 
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V. Practical Considerations and Simulations 

A. Efficient and Robust Feedback Protocol 

Mobile user equipments usually have limited computing capabilities, therefore, most systems 
require that the complexity at the user side is as low as possible. Hence, solving the full rate 
approximation problem (i.e. the minimax problem (fT3] )) may not be feasible. Fortunately, our 
analysis in Section |IV] yielded another distance function which can be used at the user side to 
uniformly bound the rate approximation error Ai?RA (fT2l) and is given by 



ds(h,i/) 



max \\{h, w)[ 

■weC 



iy,w) 



(19) 
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Fig. 4. Schematic comparison of the GDI computation at the user side; for the efficient RA distance (bottom) and chordal 
distance (top). 

Figure |4] (bottom) shows a flow chart for the CDI computation under the proposed feedback 
protocol using the new distance function. The terms |(i/j,iyj)|, for i = 1, . . . , |C| and j = 
1, . . . , |V|, only need to be computed ones and can be stored in the memory of each user. 
Hence, during the feedback phase user m must only compute \{hm, Wi)\ for all z = 1, . . . , \C\ 
and the difference 1 1 {hm, Wi) \ — \ {uj, wi) \ | for all i, j. Figure |4] (top) shows the steps that need 
to be performed to compute the CDI based on the chordal distance. To compute the chordal 
distance each user must compute \{hm, Wi)] for alH = 1, . . . , |C|. We observe that determining 
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the CDI under the proposed efficient feedback protocol is slightly more complex than using the 
minimum chordal distance. But, as we show next the proposed protocol achieves a huge sum 
rate gain. 



B. Simulations 

In the simulations we consider 3 base stations located in 3 adjacent cells and 30 users uniformly 
distributed over the network area, i.e. a radius of 250 meter around the center of the base stations. 
The physical layer is configured according to LTE [|T5l . The base station are equipped with n^ = 4 
transmit antennas and each user is equipped with ra^ = 1 or ra,. = 2 receive antenna (specified 
in the caption). The transmit codebook and feedback codebook is given by the LTE codebook 
defined in ifTSl which has A^ = 16 elements and, hence, we require i? = 4 bit to feedback back 
the CDI. The channels are modeled by the spatial channel model extended (SCME) |[T6ll using 
the urban macro scenario. The simulation parameters a summarized in Table HI 



TABLE I 

Simulation Parameters 



Parameter 


Value/Assumption 


Number of base stations 


3 


Frequency reuse 


full 


Number of users \U\ 


30 (uniformly distributed) 


Number of transmit antennas nt 


4 (uncorrelated) 


Number of receive antennas Ur 


1 or 2 (uncorrelated) 


Receiver type 


maximum ratio combining 


Maximum number of scheduled users 


4 


per scheduling block Us 




Equivalent SNR 


153 dB 


LTE carrier frequency / bandwidth 


2 GHz / 10 MHz 


Number of PRB 


50 


Scheduling block size 


1 PRB = 12 subcarrier 


LTE channel model 


SCME (urban macro) 


Inter cell interference modeling 


explicit 



Each user reports a feedback message to the base station with maximal receive SNR. We will 
consider four different feedback strategies: 
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1) Perfect (average) CSIT: the base station knows the channel averaged over all subcarriers 
perfectly, i.e. 

Mm = — / ^ -timjy 
/=1 

where Hmj is the channel of user m on subcarrier / and F is the number of subcarriers 
for which feedback is generated. 

2) Minimum chordal distance: user m determines its CDI feedback by minimizing the chordal 
distance © to the channel hm = {u,Hm), where u is chosen to maximize \{u,Hm)\. 

3) Rate Approximation as described in Section [ni] with the rates 

F ^ 

4) Efficient Rate Approximation as described in Section IV-Al where h^ is given by the 
average channel hm. as defined for minimum chordal distance above. 

The base stations run a local and independent scheduler. In each transmission interval the 
scheduler can select Ug = 2 users for transmission on the same time and frequency resources. 
The scheduling is performed in a greedy fashion according to dBJ. For simplicity we assume no 
delay in the CSI report, scheduling, transmission or performance evaluation. 

Figure [5] depicts the CDF of the spectral efficiency for users with n^ = 1 receive antennas. The 
ZF scheme is implemented according to JHl. The PU2RC scheme is based on the same transmit 
codebook as RA and is implemented according to ifTTl . We observe that with perfect CSIT 
ZF outperforms greedy scheduling with a fixed codebook. With partial CSIT the RA scheme 
significantly outperforms ZF with a gain of approximately up to 70%. Remarkable is also the 
gain of about 35% of RA over PU2RC. Moreover, Figure [5] shows that RA with the efficient 
distance function (fT9l ) performs very close to the full RA scheme. 

Figure [6] depicts the CDF of the spectral efficiency for users with rir = 2 receive antennas. 
We observe that with perfect CSIT ZF outperforms greedy scheduling with a fixed codebook. 
With partial CSIT the RA scheme significantly outperforms all other schemes and achieves a 
gain of approximately 50% over ZF. Remarkable is also the 35% gain of RA over PU2RC. 

VI. Conclusion 

We analyzed the performance of RA and showed that it outperforms ZF for a large fraction of a 
practically relevant SNR range. Hence, a remarkable result is that it is often much better to reduce 
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ZF; partial CSIT; (mean spec. eff. 1 1 .42 bits/s/Hz) 

RA simple; partial CSIT; (mean spec. eff. 18.77 bits/s/Hz) 

RA; partial CSIT; (19.15 bits/s/Hz) 

ZF; ideal CSIT; (mean spec. eff. 33.50 bits/s/Hz) 

UB; ideal CSIT; (mean spec. eff. 27.29 bits/s/Hz) 

PU2RC; partial CSIT; (mean spec. eff. 14.15 bits/s/Hz) 



20 25 30 

net spectral efficiency [bits/sec/Hz] 



50 



Fig. 5. CDF of spectral efficiency; Comparing PU2RC and ZF under perfect and partial CSIT; Setup: rir — 1. 

flexibility at the base station in favor of having more reliable CSIT. Several advantages of RA not 
been accounted for in this paper and might be part of future work. The RA scheme separates 
the feedback and transmit codebook, which allows to design both codebooks under different 
constraints. For example, the transmit codebook could be designed to minimize the inter cell 
interference where the feedback codebook could be designed to minimize the rate approximation 
error. Furthermore, each user could load different feedback codebooks if its environment changes. 
Finally, having approximations of the individual user rates available at the base station can be 
beneficial in multi cell systems with cooperating base stations. 
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■ + - ZF; partial CSIT; (mean spec. eff. 21 .42 bits/s/Hz) 

■ © - RA; partial CSIT; fixed CB; (mean spec. eff. 31 .67 bits/s/Hz) 
■ ZF; ideal CSIT; (mean spec. eff. 46.50 bits/s/Hz) 

-©— UB; ideal CSIT; (mean spec. eff. 38.89 bits/s/Hz) 

■ B - PU2RC; partial CSIT; (mean spec. eff. 23.51 bits/s/Hz) 
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Fig. 6. CDF of spectral efficiency; Comparing RA, PU2RC and ZF under perfect and partial CSIT; Setup: rir 

Appendix A 
Proof of Theorem [H 

Proof: For any mapping of beamforming vectors tx we have 



Ai?csiT < ^H [RzF (H) - R (tt, H)] 



E 



H 



Y, log (1 + \l<Pt,ZF) - log 1 + 



\2 J^h 

mr'm,-K(m) 



m£S 



1 + A2 1 



E 



H 



m,7r(m.) i 
^m.,7r(/)J 



E 



H 



'"" V " ^^^ 

^ , A , ^2 (1 + ^-) <ZF - (1 + Kn.^t,ZF) ^1m^) 

2^ log 1 + A^ 



m£S 



l + \l 



Note that under UB tt : 5 — )■ C is simply a permutation on the beamforming vectors in C 
since U = S = {1,2, ... ,nt}. Since, 0^ ^. ^ and (p^ ^.f ^^ independent random variables from 
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Jensen's inequality we get 

Eh Y, loj 




1 + A2 



A^) Eh [0:^,zf] - (1 + A^Eh [0f^,zF] ) Eff 



Lh 



1 + A^ 

Due to the ZF criterion, for all i, h^ and lu^i^j are independent isotropic vectors, see e.g. 
(si. Hence, we have that (p'^zF — \{hrn,t^zF,rn)\'^ is a beta distributed random variable with 
parameters 1 and Ut — 2 and expectation 

The expectation Ei:/[0^^. J = Ejf [|(h,m, ''^^(m))^] is harder to obtain; we resort to bounding 
it. Suppose that without loss of generality /ii > /i2 > ••• > fi-at- Since, hm is isotropically 
distributed on §"*"^ we can chose C as the columns the standard ONB in C"*. Now a reasonable 
(but clearly suboptimal) choice of n can be obtained by processing the users in a greedy fashion, 
i.e. user m gets the beamformer 

7r(m) = aigmax\{hm,Wj)f , 

where we defined the set of possible beamformer indeces as J{m) = iS\{7r(l), vr(2), . . . , 7r(m — 
1)}. Consider now a fixed m. Then, in order to get the image of ir (m) we will have nt — m + 1 
degrees of freedom due to our assumption of non-random channel gains. 

In order to get an explicit expression for EhU {hm, ■w^7r(m)) |^] we will require the order statistics 
of elements of vectors uniformly distributed on the unit sphere C"* . More specifically we need 
expectation of iaaxj^j\hm.j\'^, where hm G C" and J' C {1,2, ... ,nt} with k = \J'\ < rit. 
From standard results we have the inequalities 



Pr ( max \hr^J' < ^/ {1 - e) log {\ J {m)\)] 

Prf||^™||2>(l + e)v^) < — 



) 



ertt 
where e > 0, so that we arrive at the following expression for the expectation 

[1 - e)\og{nt ~ m ^ \) ( ^ ^_(„,_^+i)^ _ ^ 

which proves the claim. D 



(1 + e} Ut \ 
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Appendix B 
Proof of Lemma [H 

Proof: Using W = 5 = {1, 2, . . . , ut} we get 



rm{ll, Xmhm) - ?^m(7r, ^ml^m) = ^Og 



,(^^ + i)(i + a^.(i-<.m)). 



= log = . = log IH '-^ — '—^ 

\ 1 - X„^(^)'^ , J \ 1 - A™0'* , , 

Here, we have set 

^^ ~ \^ 






^2^1' '" A2 + 1 

Similarly, the negative term can be bounded from above by 

, / "'T'rn. Trim.) | ^ , / i , '"T'm,7r(r?i) ^''■^'m,Tr{m.) 
— log = ^^ < log 1 H ^^^ -^ 

\ 1 - A,„0'^ , J V 1 - ^^^^'^ , , 

\ ^''■T^m,T7{m) / \ "'■T^ m,TT(m) 

Hence, we have 

/ A^0^ , ^ - ^rn(l)'' f A " 

/ \ u \ /Q \ ^ Till '"^m,7r(m) rn^m,TT{m) i /on\ 

'"mlTT, A™/i.;„) - r„,(7r, fi'mi^m) < max log 1 H ^ — | (20) 

max log 1 H -^ -^ (21) 

1 - 7? (/)'" 



for any i^ G V. Assume the following strategy: to minimize (1201) take Vm G V such that, if 
\{hm,w*)\'^ is maximized for w* , then Ki^m, tf;*)p > \{hm,'w*)\'^. Since C C V such vector 
always exists. By this strategy we can set 

l + A^(l-|(/i,„,^*)|2) 1 + ^^^(1-|(^„,^*)|2)' 
and we get after some calculations 

~ A^, \{h^,w*)\^ _~ \{h^M)\^ 

and finally 

^ max log I i -|- Ara I = — ~ ~~ + 



max log 1 H -^ —^ = max log 1 + 
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1 I 1 , "'T^mjTrlm) ^^^in,Tr{m) i , i i -i 

max loff 1 H ^-i —'- < max log I 1 
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~^Jj^^\{''n.,w)\'-\j{h^,w) 



< max log I 1 + A^ 



{hra,w)\^ - \{Um,w)\^\ 
1 - X„MUm.,w)\'^ 



{Um,w) \\{U,W*)\'^ 



{hn^,W*)\'\ 



{U„,,W*) l-~X„,\{Um,w)\^ 



and since \{hm,w*)\'^ < \{u,w*)\'^ we have \{um,'w)\'^ < 1 Viu E C,w ^ w* which proves the 
claim. D 



Appendix C 
Proof of Lemma [2] 

Proof: According to (fT4l) . RA aims on minimizing max^ |?"m(7r, X„ih„i,) — r.m (tt, -i^i/) | over 
the elements of V. Therefore, from (ISO)) and (|2T1) and Jensen's inequality we have 



nt 



Ai?RA < 5^ log 1 + Ejj 



m,=l 



mm max 

l>t5m>0 TT 






(m) 



+ J] log 1 + Ejf 



m^S 



mm max 

l>t5m>0 TT 



1 - X^ch'^ , . 



1 - ^m(p'' , ^ 
ii^y m,iT{m,) 



(22) 



Let us re- write the first term on the right side of (l22l) . We first exploit that whenever maxTr |</'J^ ^z^-, | > 
1 — e, e < 0.5, then by Lemma [H the error can be uniformly bounded from above by 



Xm^ 



Xle Xle 

< 



1-A^e l + A2,(l-e)- l + A^e' 

and since clearly max^ Ic^fj, , J > — and (1 — e) < —^ for e < 1 ^ we have for 

max^ |0j^,,(„)| > max (0, 1 - e) 

AmC ^ A^e ^ A^e 

1 - A^e ~ 1 + A^,;^ X + A,^ 



rit — l 



for any e > (even that for e > 1). On the other hand, we have for max^r |<^^^(^J < 

max (0, 1 — e) 



Am0m^^(m) '^rn4>m,TT{m) ^"i^m,TT{m.) ^™*A'm,7r 






< 



< 



T(m) 



1 ~ Am + A^e 



1 + A„ 



nt — l 
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Hence, we can write for some pair hm, i^ 



^^^m.,n{m) 



< 



max <; A^e, 



>m.Y'rrt,7r('m) ^i^T^m,iT 



(m) 



1 + A,- 



rit-l 



and setting 



AmC 



mm max 

l>tj,n>0 T 

vev 



'^m0m,iT{m) '"rn(Pm,'iT{m) ~ '^■rn'Prn,n*{m) "m'Pm,-K*{m) 



yields 



mm max 

l>i5m>0 T 






(m) 



1 Am0^^^(^) 



< 



Am0m,7r*{m) ^m0m,7r* 



(m) 



1-A. 



nt~l 



'^^'rm,TT*{m) '^m'Pm,-K* 



(m) 



Equivalently, for the second term on the right side of (1221) we have 









< 



max <; A^e, 



Am0m,7r(m) '^m(Pm,n{m) j 






(m.) 



max < AmC 



< 



Am(^'' ^ ^ - ^ 



(m) "^i^m,7r(m) 



l + A, 



■{■ 



AmfZ-J 



T{m) 



T{Tn) 



nt—1 



M 



Setting e = (1 + e') \Xm(pt,7T*im) ~^m0m,7r*(m)l' ^' -^ '^' ^^^'^^ *^ ^""°^ ^^"^ ^^ ^^^^^ increasing in 
e, yields 



mm max 

l>i5m>0 n 






(m) 



:i + e'i 



1-^, 



< 



Am0m,7r*(m) ''^m<Pm,7r*{m) 



^^Tm,-K{m) 



1-A. 



nt — l 






Finally, expanding the fraction with (1 — A^) and applying Jensen's inequality again proves the 
claim. D 

Appendix D 
Proof of Lemma [3] 

Proof: Without loss of generality, we assume the unitary transmit codebook C is given by 
the standard ONB and define the vectors 

0m = (</'m,7r(l)5 • • • ' ^m,-K{nt)) ' 
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that define points on the (rit — l)-simplex /C„j_i = {x E C"*~^ : x > and ||a:;||i = 1} with 
edge length \/2. In a similar manner for each element of the feedback codebook i/j G V we can 
define 

which gives N = 2^ = \V\ points Q = {q^, q2, . . . , Qj^} on the {rit — l)-simplex. Using this 
notation and i!}rn = Kn, Dra{B) can be written as 



D^ni{B) = \„^mm\\(f)'^ - q\\oo, 



which can be upperbounded by 



Dm{B) < max min \\x — q\\oo ='■ S. 



Let d := nj — 1. The idea of the proof is the following. The intersection of the balls B"^ 
with radius S and the center in JC^ is a polytope with 2d facets. The spheres 6B2 of radius 
5 are inscribed in this polytope. Next we bound the number of spheres 6B2 required to cover 
the simplex JC^K This number can be given by the covering number N{}C,5B2), which can be 
bounded from above as follows. Using the Rogers-Zong Lemma ifTSl . which states that the 
covering number N(A, B), that is, the number of convex bodies B required to cover a convex 
body A, can be upper bounded by 

N(A,B)<e(Bf-^^^^, (23) 

vol(o) 

where Q{B) > 1 is the covering density of B; if MJ^ can be tiled by translates of B then 
Q{B) = 1; if the covering has some overlap then Q{B) > 1. Further, we require the Rogers- 
Shephard inequality [[T9l . which states that 

vol(^ - i5)vol(^ nB)< r j vol(^)vol(i3). (24) 

Using the assumption vol(^n ;B) = vol(i3) we get from (|23] ) and (|24l) that the covering number 
A^(^, B) is upper bounded by 

'2d\ \o\{A) 



The volumes of the rf-simplex /C"' and the scaled £2-ball 5B2 are 

/ J _|_ 1 d/2 

vo,(^^) = ^ and vol(.Bj) = ^^^^^f, 
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where r(-) is the gamma function. Wrapping up, the covering number can be upperbounded by 



^Q^^,,f2d\Til + d/2)Vd + T 1 
Solving for 6, i.e. 



^\d d\TX<^/^ 5<^' 



2nt-2\T{l + ^)JWt 



nt~l\ /z;-\ "t-i 



B 



proves the inequality. The inequality is valid provided 5 is smaller than the inradius of the 
inscribed circle of the simplex. According to Klamkin [|20l for a regular simplex the inradius 
equals the circumradius divided by rit — l. The circumradius is easily shown by the volume ratio 



and Stirlings formula to be greater than y/rit — 1. This together with the first inequality yields 
the lower bound on B. 



U 



Appendix E 
Proof of Theorem [2] 

Proof: The terms of the sum in (fT2l) can be bounded from above as follows. 







-log 



« 



i^ -I- A2 V ^^ I ■" \ IM _|_ ,92 V rt, 
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Setting ^^ = A^ we get 






log 1 



< log 1 + 



< log 1 



l&S V^m,7r(0 r'm,n(l) 




^ , ^tA^ E;g5\{m} Ctt 



(0 






'5 1 1 + AmZlie5\{m} '/'m,7r 



(0 



riiA 



151 



+ log 1 



7 , 4'm,TT{l) ~ 4'm,TT{l) 

ies\{m} 
ntX' 



\S\ 



|5\{m}|maxU^^^ 



(m) 



»n,7r(»n) 



<21og(l + ntA^max|0^,^, 



21og 1 



Pfil 



max 



(T^ 



(m) rm,TT{m) 



m,Tr(m) r m,7r(m) 



The lower bound on — (rm(vr, Am^m) — ^m(7r,^m^'m)) can be obtained in a similar manner. 
Taking expectations and using Jensen's inequality we obtain 



E/f [rm{T^, >^mhm) " ^m(vr, A^I/„,)] < 2 log 1 + Ejf 



p^il 



a 



K 



2 ™^^ |S^m,7r(m) ~ V^m,7r(m) | 



Since max^ I'Pmnim) ~ '^m-n{m)\ dcpcuds only on the channel directions h„i it is independent of 
the channel magnitude jirn- 



Eh [rm(vr, Xmhm) - r^i-K, A„i^„)] < 2 log ( 1 
Using the RA scheme and (fT2l) yields the result. 



Prit 



E 



a^ 



H 



max 



k^ — ch^ I 

m,n{m) r m,7r(r?i) | 



D 



References 

[1] H. Weingarten, Y. Steinberg, and S. Shamai, "The Capacity Region of the Gaussian Multiple-Input Multiple-Output 

Broadcast Channel," IEEE Transactions on Information Theory, vol. 52, no. 9, pp. 3936-3964, Sep. 2006. 
[2] D. Love, R. Heath, V. N. Lau, D. Gesbert, B. Rao, and M. Andrews, "An Overview of Limited Feedback in Wireless 

Communication Systems," IEEE Journal on Selected Areas in Communications, vol. 26, no. 8, pp. 1341-1365, Oct. 2008. 
[3] N. Jindal, "MIMO Broadcast Channels With Finite-Rate Feedback," IEEE Transactions on Information Theory, vol. 52, 

no. 11, pp. 5045-5060, Nov. 2006. 
[4] K. Kiran Mukkavilli, A. Sabharwal, E. Erkip, and B. Aazhang, "On Beamforming with Finite Rate Feedback in 

Multiple-Antenna Systems," IEEE Transactions on Information Theory, vol. 49, no. 10, pp. 2562-2579, Oct. 2003. 
[5] M. Sharif and B. Hassibi, "On the Capacity of MIMO Broadcast Channels With Partial Side Information," IEEE 

Transactions on Information Theory, vol. 51, no. 2, pp. 506-522, Feb. 2005. 
[6] C. Au-Yeung and D. Love, "On the Performance of Random Vector Quantization Limited Feedback Beamforming in a 

MISO System," IEEE Transactions on Wireless Communications, vol. 6, no. 2, pp. 458^62, Feb. 2007. 



July 12, 2011 



DRAFT 



28 

[7] G. Dimic and N. Sidiropoulos, "On Downlink Beamforming with Greedy User Selection: Performance Analysis and a 

Simple New Algorithm," IEEE Transactions on Signal Processing, vol. 53, no. 10, pp. 3857-3868, Oct. 2005. 
[8] M. Trivellato, F. Boccardi, and H. Huang, "On Transceiver Design and Channel Quantization for Downlink Multiuser 

MIMO Systems with Limited Feedback," IEEE Journal on Selected Areas in Communications, vol. 26, no. 8, pp. 

1494-1504, Oct. 2008. 
[9] J. Schreck, P. Jung, G. Wunder, M. Ohm, and H.-P. Mayer, "Limited Feedback in Multiuser MIMO OFDM Systems 

Based on Rate Approximation," in IEEE Global Telecommunications Conference, Nov. 2009, pp. 1-6. 
[10] T. Yoo, N. Jindal, and A. Goldsmith, "Multi-Antenna Downlink Channels with Limited Feedback and User Selection," 

IEEE Journal on Selected Areas in Communications, vol. 25, no. 7, pp. 1478-1491, Sep. 2007. 
[11] K. Huang, J. G. Andrews, and R. Heath, "Performance of Orthogonal Beamforming for SDMA With Limited Feedback," 

IEEE Transactions on Vehicular Technology, vol. 58, no. 1, pp. 152-164, 2009. 
[12] R. de Francisco, M. Kountouris, D. T. M. Slock, and D. Gesbert, "Orthogonal Linear Beamforming in MIMO Broadcast 

Channels," in 2007 IEEE Wireless Comnumications and Networking Conference, 2007, pp. 1210-1215. 
[13] S. A. Jafar and A. J. Goldsmith, "Isotropic Fading Vector Broadcast Channels: The Scalar Upper Bound and Loss in 

Degrees of Freedom," IEEE Transactions on Information Theory, vol. 51, no. 3, pp. 848-857, 2005. 
[14] K. Boroczky, Finite packing and covering. Cambridge University Press, 2004. 
[15] 3GPP TS 36.201, "LTE Physical Layer - General Description," 2009. 
[16] D. S. Baum, J. Salo, M. Milojevic, P. Kyosti, and J. Hansen, "MATLAB implementation of the interim channel model for 

beyond-3G systems (SCME)," 2005. 
[17] R. Samsung, "DownUnk MIMO for EUTRA," 3GPP TSG RAN WGl Meeting, 2006. 
[18] C. A. Rogers and C. Zong, "Covering convex bodies by translates of convex bodies," Mathematika, vol. 44, no. 01, pp. 

215-218, Jun. 1997. 
[19] C. Rogers and G. Shephard, "The difference body of a convex body," Archiv der Mathematik, vol. 8, no. 3, pp. 220-233, 

1957. 
[20] M. S. Klamkin and G. A. Tsintsifas, "The Circumradius-Inradius Inequality for a Simplex," Mathematics Magazine, 

vol. 52, no. 1, Jan. 1979. 



July 12, 2011 DRAFT 



